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Abstract. - The nearest level spacing distribution P c (s) of d-dimensional disordered models 
(d = 1 and 2) with long-range random hopping amplitudes is investigated numerically at 
criticality. We focus on both the weak (b d 3> 1) and the strong (b d <C 1) coupling regime, 
where the parameter b~ d plays the role of the coupling constant of the model. It is found 
that P c (s) has the asymptotic form P c (s) ~ exp[— yl^s"] for s ^> 1, with the critical exponent 
a — 2 — dd/b d in the weak coupling limit and a = 1 + Cdb d in the case of strong coupling. 



It is well established that the statistical properties of spectra of disordered one-electron 
systems are closely related to the localization properties of the corresponding wavefunctions 
[1-3]. In the metallic phase, the large overlap of delocalized states, which are essentially 
structureless, induces correlations in the spectrum, leading to the well known level repulsion 
effect. If the system is invariant under rotation and under time-reversal symmetry (orthogonal 
symmetry), the normalized spacings s follow Wigner-Dyson statistics at the infinite system 
size limit: 

Pw(s) = ^sexp [-J.s 2 ] . (1) 

In contrast, in the localized regime, states with close energy levels are typically localized at 
different parts of space and have an exponentially small overlap. Their levels are therefore 
uncorrelated and the corresponding spacings are distributed according to the Poisson law 

= exp[-s] . (2) 

It has been argued that the statistics of energy levels at the disorder-induced metal- 
insulator transition (MIT) is characterized by a third universal (i.e., independent of the system 
size and of the details of the Hamiltonian model) distribution P c (s), which is different from 
both Wigner-Dyson statistics and the Poisson statistics [2,4]. The asymptotic behavior of this 
distribution for s>l has been a controversial issue and still remains unresolved. On the one 
hand, the influence of the MIT on the spectral properties was studied in Refs. [2,4] by means 
of the impurity diagram technique combined with scaling assumptions. In these studies, it 
was conjectured that 

P c (s) ~ exp[-Ks] , s > 1 , (3) 
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with k s=s 3.3, the reason for such behavior being that the Thouless energy at the transition 
point is of the order of the average level spacing (E c /A w 1), and so the levels's repulsion is 
effective only for s < 1. 

On the other hand, by mapping the energy level distribution onto the Gibbs distribution 
for a classical one-dimensional gas with a repulsive pairwise interaction, ref. [5] derived the 
following asymptotic form for P c (s): 



where the coefficient Ad depends only on the dimensionality, d, and where the critical expo- 
nent, a, which ranges in the interval f < a < 2, is related to the correlation length exponent 
v and to the dimensionality through a = 1 + (df) -1 . 

As regards the numerical description of P c (s), there is also no consensus. The exponential 
decay, eq. ©, of P c (s) has been confirmed by most groups at different MITs (see Ref. [6] 
and references therein), while an exponent am 1.2 has been found in refs. [7,8] from a fit in 
the whole range of spacings to a distribution of the form P c (s) = i?sexp[— As a ] or, indirectly, 
from the two-point correlation function of the density of states [9] . Anyway, the behavior IpTjl 
with some nontrivial f < a < 2 is what one would expect at the mobility edge. 

It should be pointed out that MITs generically take place at strong disorder (conventional 
Anderson transition, quantum Hall transition, transition in d = 2 for electrons with strong 
spin-orbit coupling, etc.). In this regime, the predicted [5] exponent a = 1 + {dv)^ 1 slightly 
deviates from unity, making it relatively difficult to see on the numerically calculated tails of 
P c (s) (e.g., at the standard Anderson transition in 3D a « 1.2). To overcome this problem, 
it is necessary to investigate transitions which occur at the opposite limit (weak coupling 
regime). This area has been left unexplored, but one would expect to find an exponent a far 
from unity and closer to the Wigner-Dyson value a = 2. 

In this work, we try to definitively solve the existing controversy about the large s asymp- 
totic form of P c (s). From results of detailed high-precision numerical investigations, we will 
show unambiguously that eq. (@J is indeed correct, while the validity of eq. © is limited to 
the case of very strong disorder (strictly at the limit of infinity coupling strength). In addition, 
we find that the exponent a in eq. Q continuously varies between 1 and 2 as the coupling 
strength of the Hamiltonian model changes from to oo. 

To this end, we performed numerical calculations of P c (s) on a generalization to d di- 
mensions of the power law random banded matrix (PRBM) model [10-19] (for closely related 
models see also Ref. [20]). The corresponding Hamiltonian, which describes non-interacting 
electrons on a disordered e?-dimensional square lattice with random long-range hopping, is 
represented by real symmetric matrices, whose entries are randomly drawn from a normal dis- 
tribution with zero mean, (Wy) = 0, and a variance which depends on the distance between 
the lattice sites n: 



We refer the reader to Ref. [19] (and references therein) for the advantages of the present 
model with respect to Hamiltonians with short-range, off-diagonal matrix elements, and for 
the many real systems of interest that can be described by Hamiltonians JSJ. 

The parameter b d in eq. JSJ is an effective bandwidth that serves as a continuous con- 
trol parameter over a whole line of criticality, i.e., for an exponent equal to d in the hop- 
ping elements Hij ~ b d [21]. Furthermore, it determines the critical dimensionless con- 
ductance in the same way as the dimensionality labels the different Anderson transitions. 



P c (s) ~ exp[-^ d s Q ] , s > 1 , 



(4) 




i = 



(5) 
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Fig. 1 - The integrated probability I c (s) of the 2D system for L = 60 at b 2 = 0.02, 0.04, 0.06, 0.08, 0.1, 
0.4, 0.8, 1 and 10 (from bottom to top). Dotted and dashed lines are 7w(s) and Ip(s), respectively. 



Each regime is characterized by its respective coupling strength, which depends on the ratio 
((|7iii | 2 }/(|7ijj j 2 )) 1 / 2 oc b~ d between diagonal disorder and the off-diagonal transition matrix 
elements of the Hamiltonian [22]. 

We remind the reader that in the two limiting cases of the ID model, 6^1 and b 
1, which correspond to the weak- and the strong-disorder limits, respectively, some critical 
properties (spectral compressibility, correlation dimension, ...) have been derived analytically 
by mapping Hamiltonian (JSJ) onto an effective u-model of a one-dimensional nature [10-14]. 
We stress that, unlike the ID PRBM model, it has not until now been possible to analytically 
solve the 2D disordered models with long-range transfer terms. 

The system size ranges between L = 1000 and 4000 in ID, and between 20 and 100 in 
2D, whereas b d ranges the interval 0.02 < b d < 10. We consider a small energy window, 
containing about 10% of the states around the center of the spectral band. The number of 
random realizations is such that the number of critical levels included for each L is roughly 
1.2 x 10 6 , except for the larger system size in 2D, for which this number is about 3 x 10 5 . In 
order to reduce edge effects, periodic boundary conditions are included. 

For the computation of P c (s), we unfold the spectrum in each case to a constant density, 
and rescale it so as to have the mean spacing equal to unity. In order to diminish the magni- 
tude of the relative fluctuations and to analyze the asymptotic behavior in detail, it is more 
convenient to consider the cumulative level spacing distribution function J(s) = J P(s')ds' . 
Note that the integration does not change the asymptotic behavior of P(s). The Wigner 
surmise, eq. JPl, and the Poisson distribution, eq. @, yield Iw(s) = exp[— 7rs 2 /4] and 
Ip(s) = exp[— s], respectively. 

Figure^displays our results for the integrated probability I c (s) of the 2D system for L = 60 
at b 2 = 0.02, 0.04, 0.06, 0.08, 0.1, 0.4, 0.8, 1 and 10, which are depicted consecutively from 
bottom to top. Dotted and dashed lines, which correspond to Iw(s) and Ip(s), respectively, 
are given for comparison. A gradual crossover in the large s tail of I c {s) from the Poisson 
to the Wigner-Dyson limiting forms as one increases the inverse coupling constant b 2 of the 
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Fig. 2 - Log-log plot of the integrated probability I c (s) of the ID system at b = 0.02, 0.1, 0.2, 0.5, 
1.0, 1.5, 2 and 5 (from bottom to top) and different system sizes L = 1000 (circles), 2000 (squares) 
and 4000 (diamonds). Dotted and dashed lines are Iw(s) and Ip(s), respectively, and the straight 
lines are fits to eq. @. 

Fig. 3 - The & d -dependence of the critical exponent a for the ID (circles) and 2D (squares) disordered 
systems. Solid lines are fits to eqs. © corresponding to the limiting cases of weak (b d S> 1) and 
strong (b d <C 1) disorder. 



model can clearly be seen. So, we can therefore expect an exponent a in eq. (@J, which spans 
the interval [1,2], in agreement with ref. [5]. 

Next we consider the behavior of I c (s) with system size L. The results for s large of the 
critical I c {s) for the ID system at different values of b are shown in a log- log scale in fig. [21 for 
different system sizes: L = 1000 (circles), 2000 (squares) and 4000 (diamonds). As in the 2D 
case (see fig. one can viasualize the crossover between the small-& and large-6 asymptotics. 
Note that I c (s) is an L- independent universal scale- invariant function that interpolates, as 
previously mentioned, between Wigner and Poisson limits . This result confirms the existence 
of a critical distribution exactly at the transition. Dotted and dashed lines correspond to 
iw(s) and ip(s), respectively. We checked that the normalized variances of P c {s) are indeed 
scale- invariant at each critical point studied [23]. The straight line behavior of the data in 
such a plot at all values of b considered is undoubtedly consistent with a ^-dependent exponent 
a in eq. (QJ. The values of b reported are 0.02, 0.1, 0.2, 0.5, 1.0, 1.5, 2 and 5, from bottom 
to top. The best fit to eq. (@J in the interval 2.5 < s < 5 for small b and 2.5 < s < 4 for 
large b, yields a = 1.008, 1.040, 1.085, 1.201, 1.376, 1.600, 1.703 and 1.866, respectively, thus 
confirming the result of [5]. Note that for the large energy ranges considered, where I c (s) vary 
by one to three orders of magnitude, the quality of the fits, which are represented as solid 
straight lines, is evident. 

The disorder dependence of the critical exponent, a, as obtained from the previous fits 
for the ID (circles) and 2D (squares) systems is shown in fig. [3] in the broad range of the 
parameter b d of the PRBM model. For both dimensions, d — 1 and 2, it clearly changes 
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Fig. 4 - Log-log plot of -Pc(s) of the ID system for different sizes L = 1000 (circles), 2000 (squares) 
and 4000 (diamonds) and disorder strength b = 5, 0.1 and 0.05 (from bottom to top). Dotted and 
dashed lines are the Wigner surmise, eq. Q, and the Poisson distribution, eq. JS), respectively. Solid 
lines are fits, in the given intervals, to the form P c (s) = Cs. 

continuously from 1 as b d — > to 2 as b d — > oo. In the two limiting cases of weak (b d ^> 1) 
and strong (b d <C 1) disorder regimes it can be well fitted by 



respectively. These fits are shown as solid lines in fig. [3] The fitting parameters are a\ = 0.60, 
a,2 = 0.33, ci = 0.42 and ci = 0.76. Note that the different values of these parameters reflect 
its dependence on the dimensionality. From eq. ©, the Poissonian tail of P c (s), eq. ©, is 
recovered for large spacings at the limit of very strong coupling b d — » 0. So, we conclude that 
in the case of very strongly coupled Hamiltonians only, eq. J3J losses its validity and eq. © 
applies. 

The observed 6 d -dependence of a in eq. © is not surprising, since other critical proper- 
ties, such as the spectral compressibility or the correlation dimension di in ID, for which 
analytical treatment is feasible, present a similar behavior towards b. Specifically, \ — Xjlixb 
(6 » 1), % = 1 - 46 (6 < 1), d 2 = 1 - 1/nb (b > 1), and d 2 = 2b (b < 1) were derived in 
refs. [10, 12, 13]. We stress that eq. © is based on numerical results and at present it should 
be considered as a conjecture. So, further analytical work is needed to check this form of the 
critical exponent and its origin from the model JSJ. 

Finally, we present the limiting behavior of P c (s) as s — > 0. From general considerations 
for the orthogonal symmetry P c (s) ~ s at small s <C 1 [2,24]. The results for the ID case 
at different values of b for various system sizes, L = 1000 (circles), 2000 (squares) and 4000 
(diamonds), are plotted in fig. 4. Dotted and dashed lines are the Wigner surmise, eq. £[J, 
and the Poisson distribution, eq. respectively. Here we find that P c (s) ~ s for all disorder 
regimes in accordance with the predictions of refs. [2,24]. The slopes of straight lines fitting 
the data in the intervals shown are 1.63, 6.07, and 11.30 at b = 5, 0.1 and 0.05, respectively. 




b d > 1, 
b d < 1, 



(6) 



6 



EUROPHYSICS LETTERS 



We have checked that, for the 2D case, the same linear behavior of P c (s) towards s is fulfilled 
in the whole range considered of the parameter 6 2 . 

To summarize, we have investigated the critical level spacing distribution P c (s) of non- 
interacting electrons on a c?-dimensional disordered system with long-range transfer terms in 
the whole range of the coupling constant b~ d . P c (s) is found to be scale independent at all 
values of b~ d . The large s part of P c {s) obtained is shown to have an exp[— AdS a ] decay with 
1 < a < 2. We determined the disorder dependence of a in both the strong (b d <C 1) and the 
weak (b d 3> 1) coupling regimes. At the limit of very strong disorder b d — > 0, we found that 
a — > 1 and so we obtain the expected results of the Poissonian decay predicted in refs. [2,4]. 
The small-s behavior of P c (s) ~ s is in agreement with the analytical predictions at all values 
of b d . 
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